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Abstract. In this paper, a single sum formula for the lineariza- 
tion coefficients of the Bessel polynomials is given. In three special 
cases this formula reduces indeed to either Atia and Zeng's formula 
(Ramanujan Journal, Doi 10.1007/slll39-011-9348-4) or Berg and 
Vignat's formulas in their proof of the positivity results about these 
coefficients (Constructive Approximation, 27 (2008), 15-32). As a 
bonus, a formula reducing a sum of hypergeometric functions 3F2 
to 2F1 is obtained. 
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1. Introduction 
The Bessel polynomials g„ of degree n are defined by 

fc=0 ^ 

where we use the Pochhammer symbol {z)n '■= z{z + 1) . . . {z + n — 1) 
for 2; G C, n = 0, 1, .... The first values are 

qo{u) = l, qi{u) = l + u, g2(M) = l + M + y- 



Some recursion formulas for g„ are 

qn+i{u) = qn{u) + -^-^—^(ln-i{u), U > 1, (2) 

I'niu) = qn{u) - -qn-l{u), 71 > 1. (3) 

2n — 1 

Using hypergeometric functions, we have qn{u) = iFi{—n; —2n; 2u). 
They are normalized according to qn{0) = 1, and thus differ from the 
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monic polynomials 6'„(m) in Grosswald's monograph [1]: 




{2n)\ 



qniu). 



The polynomials On are sometimes called the reverse Bessel polynomials 
and Uniu) = u^Oni^^) the ordinary Bessel polynomials. These Bessel 
polynomials are, then, written as 



The linearization problem is the problem of finding the coefficients 
/^^."■''"■'(ai, 02) in the expansion of the product Pn{0'iu)Qm{0'2u) of two 
polynomials systems in terms of a third sequence of polynomials Rk{u), 



The polynomials -P„, Qm and Rk belong to three different polynomial 
families. In the case P = Q = R and Oi = 02 = 1, we get the (stan- 
dard) linearization or Clebsch-Gordan-type problem. If Qm{u) = 1, we 
are faced with the so-called connection problem. 

In the case P = Q = R and ai = a, 02 = 1 — a, we get the Berg- 
Vignat linearization problem. And, finally, in the case P = Q = R and 
for any 01,02, we get a new linearization problem. 

In this paper, we are interested by this new linearization problem 
and by the linearization coefficients /^^"'"^^(ai, 02) in the case of the 
Bessel polynomials which are defined by 





n+m 





k=0 



n+m 





k=0 



For example, we have 
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where 

(01,02) = 143ai^a^, 

I3j^'^\ai, 02) = — ai'a^ (12 aia2 — 5 Oi — 202) , 

/^e^'^H'^i, 0,2) = -3- cti 0.2^ (—140 ai^ 02 + 35 ai^ + 30 ai 02 + 5 02^ — 56 ai a'l + 126 ai^ 02^) , 
o 

/3f '^Vi, «2) = 02^(02^ + 42 01^02 + 28 + 84 a-i" - 84 a2 
+ 2100^02 - 21aia2 - 1470^02 + 15aia2 - 1260103), 



fif'^^ip^x^ 0'2) = - a2(245a'^a2 + 21a!^a2 — 140a^a2 — 140a'^a2 + 35aia2 — 75aia2 — ^^a\a\ 
3 

+ 210a^a2 - 210a^a2 + 5a2 + 56a^a2 + 21a5 - 5a2 + 35aia2), 

(3 5), . Q 5 Q 4 5 c 35 Q o 5 Q 50 o 75 4 20 n 
P3 (cii, 02) = ai + -a^a^ - -a^a^ - —a^a^ + -a^ - —a^a^ + —a^a^ + —a^a^ 

50/1 0/1 O OQ^O^O '^O 

02 — 10a;^a2 + 6a;ia2 + 300^^02 a^a^ + 20aja2 — 10aja2 

2 X 3 

5 c 295 

/32^'^''(ai, 02) = (fli + a2 — l)(140a?a2 + 126a? — 315a?a2 — 385aia2 + 315aia2 

105 

+ 70ai4 - 704 + 140a^ + 5a^), 

/3p'^''(ai, 02) = T7-(oi + a2 — l)(3a^ + 15aia2 — 15ai — 15a2 + 503), 

/3o^'^^(ai, 02) = -ai - a2 + 1. 
For n, m > 1 and Gx — O, 0/2 

= 1 — a. Berg and Vignat [2j have proved 
the following recurrence relation for (a, 1 — a) which they denoted 
by/3i"'"V) P Lemma 3.6]: 

2fcTT^^+^ " ^^^^ ^^^+2^;^^^ 

for /c = 0, l,...,m + 72 — 1. From ([8]) they derived the positivity of 
/3^"'™'^(a) when < a < 1 and also that /^^"'"^^(a) = for k < min(m, n). 
Recently, with J. Zeng |1], we improved this result by giving the ex- 
plicit single-sum formula for /3^"''"'' (a) which was missing in their paper 

m- 

In this paper, our main result is twofold: 

- for any 01,02, a recurrence relation for /3^"'™''(oi, 02) is given. This 
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recurrence relation reduces to the recurrence system ([H]) when ai = a 
and a2 = I — a. 

- for any ai,a2, an exphcit single sum formula for /3^"''"'''(ai, 02), which 
provides actually the unique solution of the recurrence relation and, 
then, becomes a generalization of (3^'^\o') given by Atia and Zeng in 
yy when tti = a and 02 = 1 — a. 



Lemma 1. Forn,m > 1, the recurrence relation fulfilled by I3l!^'"^\ai,a2), < 
k < n + m is given by 

and for 0<k<n + m — 1, we have 

al{2m - l)(2m + 1) ^n-i,m+i)f . 
" ai(2n-l)(2n + l) '''^ 

n(n,m^l)( x (2m - 1 ) (2m + 1) 

= ai(2n-l)(2n + l) 

Proof. In one hand we have 

n+m 

qn+i{aiu)qm~i{a2u) = I3^^^^'"'~^\ai, a2)qk{u), 

k=0 

in the other hand, using ([2]), we have, 

g„+i(aiM)g^_i(a2u) = g„(ain) + — -\— — — -g„_i(ain) gm-i(a2^^) 

V (2?2 - l)(2n + 1) y 

= g„(aiM)gm_i(a2M) + -— -^-7; — — — g„_i(aiM)gm_i(a2n) 

[2n — l)[2n + 1) 

, , , , af(2m - l)(2m + 1) , , a^u"^ , . 

= g.(a,n)g_,(a2«) + ,2(2^ _ ^ 1) gn-i(^i^) _ ^^(^^ ^ 1) ^^-1(^2^) 

a^(2m — l)(2m + 1) / 

= qn[a'\u)qm~\{a2u) + ^277^ — — — — T-g„-i(aiM) gm+i(a2«) - qmia2u) 

a2{2n — l)(2n +1) \ 

where we used again (j2]), finally, we obtain 

a?(2m - l)(2m + 1 

ai(2n- l)(2n+ 1) 

a2(2m-l)(2m + l) . n / n 
= '?n(aiM)gm_i(a2M) - ^2^-2^ _ i'j(2n + 1) '^"-^('^^^^^'"("^w), 

and because of the degree of polynomials qk{u) we have ([9]) and for 
0<A;<?7, + m — Iwe have ffTOl) . 



g„+i(aiM)gm_i(a2M) - ^2/o„ , i\ '?n-i(Qiw)gm+i(a2w) 
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Theorem 2. For i = 0, 1, n + m, we have 

n{n,m) / \ 



—m+k — n+fc/1 irt\ 

02 (l/2)fc 



4m+n-fc(^ + n - A;)!(l/2)„(l/2)„ 

m+n—k 

J2 a^+""*^"^(r^""'')(n + 1 - i)2i 

i=0 
m+n—k—i 

j=0 

(11) 

which we write using 3F2 hypergeometric functions as 
Theorem 3. For i = 0,1, ...,n + m, we have 

n{n,m) / \ _ U-^ U2 l^-L/^JA; 

^'''''''^ " 4-+-^(m + n-fc)!(l/2)„(l/2)^ 

m+ra— fe 

5^ 0^l{ZXT-k-^)^-^ + k + l + l)2im+n-k^i){m - I 
1=0 

, k + 2, —k — 1, —i . i 

Remarks. 

1. This formula was deduced using the same approach done in [1] pages 
4 and 5 by, just, changing a by ai and 1 — a by 02- 

2. To compute this formula with, for example. Maple, one should com- 
pute /^i'^^^ (01,02), /3^'^^^_i(ai,a2),..., f3^''''^\ai,a2) and then replace 
n, m by their values (please see the Maple program given in the end of 
this paper). 



Proof of theorem 3. Let us, first, prove that (IT^ fulfils 

n{n+l,m—l) / \ 
Pn+m 1^1, O2 J 



(n+i,m-i), „ , _ a^iar'v^r(l/2 + n + m) 



r(n + 3/2)r(m- 1/2) 
and 

, , _ aria^+^ V?r(l/2 + n + m) 

Pn+m K'^l) 0,2) — 



r(n- l/2)r(m + 3/2) ' 
then 

a^(2m-l)(2m + l) ^ a^^a^ 'v^r(l/2 + n + m) 

ai(2n - l)(2n + 1) ^"^"^ ""'^ r(n + 3/2)r(m - 1/2) 
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Second, we prove that ( fT2l) fulfils ( fTOjl . so let us substract the rhs 
from Ihs of f lTU]) to obtain 

o{n+i,m-i), X ai(2m - l)(2m + 1) 

~ ai(2n-l)(2n + l) ''^^ 

^(n,m-i). a^(2m-l)(2m+l) 
+ ai(2n-l)(2n + l) 



i^"+i+'=a^-^v^r(l/2 + A;) 



■X 



4™+"-^'(m + n - A;)!r(n + 3/2)r(m - 1/2) 

(m+n—k 
«lC"=t.)(-^ + A; + ^ + 2)2(„+„_fe_.)(m - z)2. 
i=0 

3F2([-i, k + 2, -k - 1], [-m + 1 - z, m - z], 02)0^* 

rrt+n— fc 
i=0 



iF2{[—i, k + 2, —k — 1], [—m — i — l,m — z + 2], 02)02 



-m+l+fc m— 1 



-v^r(i/2 + fc) 



4m+n-fe-i(^ + n - A; - l)!r(ra + l/2)r(m - 1/2) 

(m+n—k— I 
Yl <Cln-l-ll)i-^ + k + Z + 2)2in+m-k-^-l){m - 2)2. 
i=0 

3F2{[-i, k + 2,-k - 1], [-m - i + l,m- z], 02)0^* 

( _|_ 1 /0\ m +n-k -l 



3F2([-2, k + 2,-k - 1], [-m - i,m - i + 1], 02)02 ^ , 

because 

(2m-l)(2m + l) 1 



(2n - l)(2n + l)r(n - l/2)r(m + 3/2) T{n + 3/2)r(m - 1/2) ' 
Cancelling the common factor 

+ k) 

4^+n-k-i(^^ + n - A; - l)!r(n + l/2)r(m - 1/2) 
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in both quantities, we get 

1 



4:{m + n- k){n + 1/2) 

/■m+n—k 

( «lGn:tJ(-"^ + ^ + ^ + 2)2(n+m-fe-0(m-i)2i 

^ i=0 

3-p2([-«, k + 2,—k — 1], [—m + 1 - i,m - i], 02)02' 

m+n— fc 

- ^ aiC+;;:ti)(-m + /c + ^)2(n+rr^-/c-^)(m-^ + 2)2, 
1=0 

3F2([-i, /c + 2, -/c - 1], [-m - i - l,m - i + 2], 02)03 

^m+n— fe— 1 

- ( «i C^Itnlfcl'-i) (-^ + A; + i + 2)2(„+^-fc-i-i) (m - i)2i 

3F2([-i, k + 2,—k — l\, [—m — i + l,m — i], 02)03 * 

" (11 1/2)02 "^^^ «lCn:^='-l)(-^ + ^ + ^ + l)2(n+rn-.-.-l)(m - ^ + 1)2. 

3F2([-i, A; + 2, -A; - 1], [-m - i,m - i + 1], 02)03'^ , 
equivalently 

1 

■X 



4{m + n- k){n + l/2) 

(m+n—k 
<CVn-U)i-^ + k + i + 2)2in+m-k-^){m - t)2i 
i=0 

3F2([-i, A; + 2, -A; - 1], [—m + 1 — i,m — i], 02)03* 

m+n~k 

- <CTn-Ui-^ + k + ^)2(n+rn-k-^){nl-^ + 2)2^ 

i=0 

3F2([-i, A; + 2, -A; - 1], [-m - i - l,m - i + 2], 02)03 



m+n— fc— 1 

Y + k + i + 2hin+m-k-^-l){m - th 



3F2{[—i, A; + 2, -A; - 1], [—m — i + l,m — i], 02)0. 



I 

2 
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{m + l/2)ai 

z 

=0 

3F2([-i, /c + 2, -/c - 1], [-m - i.m - i + 1], 02)03 



To prove that this expression vanishes, it suffices to prove that the 
coefficient of a{ vanishes. The coefficient of a{ is given by 



■X 



A(m + n- k)(n + l/2) 
(^a{{ZX:-'k-jK-m + k + j+ 2)2in+m-k-j){^ - 3)23 

3F2{[-j:k + 2,-k - l],[-m + l - j, m - j], 02)02^ 
-a{{ZXT-l-j)i-m + k + jhin+m-k-dm - j + 2)2,- 

3-?^2([-J, k + 2,-k-l], [-m - j - 1, m - j + 2], 02)02"- 



3-^2([-i, A; + 2, -A; - 1], [-m - j + l,m- j], 02)0^-^ 

+ ^^^^«r^C:n:^:])(-^ + k + Jhin+m-k-d^ " J + 2) 

3F2([-0- - 1), A; + 2, -A; - 1], [-m - (j - 1), m - (j - 1) + 1], 02)0" ^^"'^ 



2(i-l) 



A short computation (with Maple) of this quantity gives zero: 

Ql := ((ol* * hinomial{n + m — k,n + m — k—i) * pochhammer{—m + 
A: + i + 2,2*ri + 2*m — 2*A; — 2*i)* pochhammer{m — i,2 * i) * 
hyper geom{[—i, k + 2, —k — 1], [— m — i + l,m — i], a2) * a2(~*^) — (ol* * 
binomial{n + m — k,n + m — k — i)* pochhammer{—m + A; + i, 2*n + 
2*m — 2*A; — 2*i) *pochhammer{m — i + 2, 2*i)* hyper geom{[—i, k + 
2,-k-l], [m-i + 2, -m-i-l],a2)*a2(-*)))-4*(n+m-A;)*(n + l/2)* 
(aP * binomial{n + m — k — l,n + m — k — i — 1)* pochhammer{—m + 
A; + i + 2,2*n + 2*m — 2*A; — 2*i — 2)* pochhammer{m — i,2*i) * 
hyper geom{[—i, A; + 2, —k — 1], [— m — i + 1, m — i], a2) * o2^~'^); Q2 := 
— 4*(n+m— A;)*(m+l/2)*al/o2*(al''*~"^-'*6momia/(n+m— A;— 1, n+m— 
k—{i—l) — l)*pochhammer{—m+k+{i — l) + l, 2*n+2*m—2*k—2*{i — 
l)—2)*pochhammer{m—{i — l) + l, 2*{i — l))*hypergeom{[—{i — l), k + 
2, -k-1], [-m-(i-l), m-(i-l)+l], o2)*o2(-^+i)); simplify{Ql-Q2)- 
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Particular case. 

Let us prove that ffTUj) reduces to when ai = a, a2 = I — a. 
From ( !T0|) we have 

PI '(aA — a) — —pi Ma, 1 — a 

i'^ - n{n+l,m-l) , ^ (2m + 1 ) N /i Q^ 

= 2m - 1 (4n2-l) ^(a, 1 - a). (13) 

equivalent ly 

a^T^^' '"-^""'^ (4tf-i) 

Adding ■:2^zi(^k'~^'^\^i 1 — o) to both sides, we get 

(1 n{n,'m.—l) / \ , n(n—l,m)f \ 

fa, 1 — a) 

2^ - 1) ' " ^ ' ' 2m - 1 ' " 

According to ([8]) le Ihs is equal to 2fc+T'^i+r''('^' I — a). 
Using (Q, the rhs becomes 

a^(2m+l) , , , , ,^ 2m + l, , , , , 

- (4^2 _ (ln~i{au}qm+i{{l - a)u} + ^ (1 + 2^ + 1 )Qn~i{au}qm{{l - a}u} 

n '^^^^^^^^^^^^^^^-^ _ ^^^y 

2m — 1 



a^(2m+2l«("-l''»)^^ 1 _ ^ ~ n(n+\,m~\). 

(4^2-1)^'= ^""'^ ''^+2m-l^'= * 



Using ([2]), we obtain 



a2(2m+l) I ^( ((. ^ ^ , (l-a)V ..^ . . 



(4n2-l) ""^^^ 4m2-l 



a^ 2m +1, , , , ,^ , , 

^ 2n-\ ^ ^ 2n + 1 )g"-H^^)g"^((l " ^)") 



(l-a)V aV 



2m- 1 V 4n2 - 1 

After simplification, we get the rhs of (IH 



qn{au) + — ^ -gn-i(aM) gm-i((l - a)^)- 
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2. Applications 

1- These coefficients /3j^''^\ai, 02) with ai + 02 7^ 1 can be apphed in: 
if X and Y are two student random variables with n and m degrees of 
freedom then the hnear combination aiX + has for characteristic 
function 

n+m 

e(-^--^-)q^{a,u)qUa2u) = e^-^^-^^^^ /^^'"(ai, a2)gfc(n), 

fe=0 

On the other hand, we have 

aiX+a2Y = (ai + a2)i—^ — -X+—^ — -Y) = (ai + a2)idiX+d2Y) 

(01 + 02) (01 + 02) 

with o 1 + 02 = 1 then it exists a NON TRIVIAL relation between the 
coefficients /3^"'™^(oi, 02) and the coefficients /3(oi, 1 — oi) which is not 
clear in their expressions. 

2- For oi = o, 02 = 1 — o, these coefficients /3^'"^\a,l — a) give a 
formula reducing a sum of hypergeometric functions 3-F2 to 2-^1: 

Theorem 4. Taking into account f lT2|) and formulas (7) — (8) given in 
|1] , we get: for k > [(n + m — l)/2] 

^2«+2m-2fc(^ - o)'"-"+^r(n + m + 2) m + A; + l, -2n - 2m + 2A; . 1 . 

r(-n - m + 2A; + 2) ^ - m + 2A; + 2 'a 

m+n— 

= Yl AZX:-l-^)i-rn + k + ^ + l)2irn+n-k-^){rn-^ + lh 



. , , - J2i 

+ 2, -k-1, -i 
-m — i, m — i + 



3F2( ' ' . ;l-o)(l-o) 



and /or k < [{n + m — l)/2j 

(-a)"+^+'"(l - o)-"'-"+^r(2n + 2m - 2A: + l)r(n - k) ^ n - k, -n - m - 1 _ I, 

m+n—k 

k + 2, -k-1, -i 

— m — z, m — i + 



3^2( 1 ^' , ^ ;l-o)(l-o)- 
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Please find next a Maple program which, not only, tests that our for- 
mula is right from min{n, m) to n + m but, also, show that 

f3^j^'"^{a, 1 — a) = for /c < ■min{n, m). 

> restart; 

> A := (n, m)— > q{n, al*u)* q{m, a2*u) — sum{beta{n, m, k, al, a2) 
* q{k, u),k — min{n, m)..n + m :) 

We assume n less or equal m. This program runs from min{n, m) untill 
n + m, take any values of n, m, for example 2 and 8 

> AA := A{2,8) :; 

> alpha ■.= {n,k)->n\*{2*n-k)\* 2^1(2 * n)!/(n - k)\lk\ : 

> q :— {n, u)— > sum{alpha{n, k) * u'', k — 0..n) :; 

> beta := (n, m, k, al, a2)- > /actor (al * a2™ * Pi^^/^) 

*GAMMA{l/2 + k) * sum{aV * hinomial{n + m — k,n + m — k — i) 
*pochhammer{—m + k + i + l,2*n + 2*m — 2*k — 2*i) 
*pochhammer{m — i + 1,2 * i) * simplif y {hyper geomt{n, m, i, k)) 
*a2("*\i = 0..n + m - k)/{A''n + m - k))/{n -\- m - k)\ 
/GAMMA{n + l/2)/GAMMA{m + 1/2)) :; 

> AAA := f actor {AA) :; 

> hypergeomt := {n, m, i,k)— > simplify {hyper geom{[—i, k + 2, —k — 1], 
[—m — i,m — i + 1], a2)) : 

> collect{f actor {simplif y {AAA)) , u); 

1/5 * (—1 + al + a2) * (5 * a2 * al — 5 * al — 5 * a2 + 2 * a2^) * u 
+1/5 * (-1 + al + a2) * (5 * a2 * al - 5 * al - 5 * a2 + 2 * a2^ - 5); 
We meet again that beta{n, m, k,a,l — a) vanish for k < min{n, m). 



